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We study the role of transverse spatial degrees of freedom in the dynamics of signal-idler phase 
locked states in type-II Optical Parametric Oscillators. Phase locking stems from signal-idler po- 
larization coupling which arises if the cavity birefringence and/or dichroism is not matched to the 
nonlinear crystal birefringence. Spontaneous Bloch domain wall formation is theoretically predicted 
and numerically studied. Bloch walls connect, by means of a polarization transformation, homoge- 
neous regions of self-phase locked solutions. The parameter range for their existence is analytically 
found. The polarization properties and the dynamics of walls in one- and two transverse spatial 
dimensions is explained. Transition from Bloch to Ising walls is characterized, the control parameter 
being the linear coupling strength. Wall dynamics governs spatiotemporal dynamical states of the 
system, which include transient curvature driven domain growth, persistent dynamics dominated by 
spiraling defects for Bloch walls, and labyrinthine pattern formation for Ising walls. 
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I. INTRODUCTION 



Optical Parametric Oscillators (OPO) are versatile nonlinear optical devices j^] with a variety of possible applica- 
tions including useful alternatives to lasers and the generation of light with non classical properties ||,^. For optical 
cavities with large Fresnel number, they have also become a paradigm for the study of transverse pattern formation 
that arise in optical systems as a consequence of diffraction and nonlinearity |4| |^]. Experimental observations of 
such patterns have been reported 0. Recent interest in these transverse structures in OPOs arise from the study of 
macroscopic manifestations of quantum phenomena in the spatial correlations present in these patterns |^ , as well 
as from the study of spatially localized structures, such as domain walls and cavity solitons fsHlSl], with possible 
applications in all-optical signal processing. 

In type-I OPO the signal and idler fields generated in the down conversion process have the same state of linear 
polarization. In type-II OPO signal and idler are orthogonally polarized. The additional vectorial degree of freedom 
of type-II OPO is very interesting from the point of view of possible new nonlinear phenomena. An interesting 
example of these new possibilities has been recently observed experimentally and described theoretically ]l6| , p^ when 
considering a direct intra-cavity polarization coupling: It is possible to reach a situation of frequency degeneracy 
and phase locking between between the orthogonally polarized signal and idler fields. This is important because, 
without the direct polarization coupling, a type-II OPO remains non degenerate at frequency degeneracy because of 
polarization. In this phase locked situation the polarization of the output field is determined by the locked value of 
the relative phase between signal and idler, which can be tuned by changes of experimentally accessible parameters. 
This device has been proposed as a candidate to generate bright quantum entangled states. Our general aim in this 
paper is to consider such phase locked states in a cavity of large Fresnel number, and to explore how the transverse 
spatial degrees freedom enter in the description of the phenomenon. We find that equivalent phase locked solutions 
grow locally, forming spatial domains separated by domain walls. We study the nature and dynamics of these domain 
walls. 

When considering transverse spatial degrees of freedom in a type-II OPO without direct polarization coupling there 
are two different regimes. In one of them, characterized by an effective negative detuning [ p^ , a finite wave number 
is selected at threshold. In a second regime to be considered in this paper, and which occurs for the opposite sign of 
detuning, homogeneous solutions are selected at threshold. However, there is a continuous of possible solutions with 
arbitrary relative phase between signal and idler. Therefore, there are no possible domain walls. This is different 
of what happens in type-I OPO in which, for the equivalent regime of detunings, homogeneous solutions with two 
possible opposite phases can be selected. As a consequence, spatial phase domains appear in the system separated 
by domain walls p|,p|,|T]] | l 5| | T9|j20|] . Such domain walls are of Ising type, that is fronts for which the field vanishes 
at the core of the wallK|7 Domain walls with the same symmetry properties have been also reported for a variety 
of other optical systems ]22|-|2^. Direct polarization coupling in type-II OPO breaks the invariance under changes 
of relative phase, allowing for the formation of domain walls. The novelty here is that these walls can be either 
of Ising or Bloch type |l0y2l| ]. Differences between Ising and Bloch walls are that there are two equivalent Bloch 
ways to connect two spatial domains (symmetry breaking) and that Bloch walls can move spontaneously, leading to 
complicated persistent dynamical states of the system. The transition from Ising to Bloch walls is controlled by the 
strength of the polarization coupling. Bloch walls have been recently predicted in other optical systems 1 29| . 

Direct polarization coupling between signal and idler in type-II OPO has been discussed in the literature | I^ , p7| , p0| , ^ 
by considering the insertion in the optical cavity of wave-plates (such as quarter- wave or half- wave) . In these previous 
studies the transverse spatial degrees of freedom were not considered. The generic phenomena described in this 
paper are expected for any form of direct polarization coupling. However, we address specifically the signal-idler 
coupling arising from birefringence and dichroism of the cavity mirrors, although our general equations give a general 
representation of possible forms of polarization coupling. A small amount of birefringence or dichroism is always 
present due to weak cavity imperfections and therefore the phenomenon considered here should be generally present 
in type-II OPO's. 

The paper is organized as follows. Section 2 presents our general model equations, which are derived in detail 
in an Appendix. In section 3 we calculate the OPO threshold, we describe the possible stationary phase-locked 
homogeneous solutions, and their polarization properties are characterized in terms of the Stokes parameters. In 
Section 4 we discuss domain walls in one transverse spatial dimension (ID): Ising and Bloch walls, their dynamics and 
the Bloch-Ising transition are characterized. We also describe the polarization properties of these walls. Sections 5 
and 6 describe the dynamical states in 2 transverse spatial dimensions (2D) governed, respectively, by Bloch or Ising 
walls. Our main conclusions are summarized in Section 7. 
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II. EQUATIONS FOR A TYPE-II OPO WITH DIRECT POLARIZATION COUPLING 



A type-II OPO that consists of a ring optical resonator, filled with a birefringent, nonlinear quadratic medium 
will be considered. The device is externally pumped by a laser beam, uniform in the plane transverse to the cavity 
longitudinal axes and of frequency Up. We take into account effects of birefringence and dichroism, that can be due 
either to small imperfections of the cavity mirrors or to weakly birefringent (e.g. wave plates) or dichroic optical 
devices inserted in the optical cavity. The derivation of the governing equations is presented, for clarity, in the 
appendix. For the sake of simplicity, the cavity birefringence and dichroism are supposed to be only due to one of the 
resonator mirrors. Note also that, in general, the mirror principal axes (i.e. those along which the Jones matrix |3^ ] 
that represents the polarization transformation is diagonal) are rotated with respect to the principal axes of the crystal 
(i.e. those along which the susceptibility matrix is diagonal). This rotation angle (0) is an important experimental 
parameter through which the strength of the effects we describe below can be controlled. We note that the equations 
are obtained in the mean field, paraxial and single longitudinal mode approximation for all the fields involved. 
The equations that describe the time evolution couple together four field envelopes that depend on the transverse 
coordinates x,y: the linear polarization components of the intra cavity field at the pump frequency: Bx^y{x,y,t); the 
signal and idler fields: Ax^y{x,y,t). The signal and the idler can be either frequency degenerate or non-degenerate, 
depending on the frequency selection rules imposed by the combined effects of the parametric down-conversion, 
the cavity resonances and the phase- matching [p3|-p5|, but they are always polarization non-degenerate (type-II 
interaction). Hereafter the freq uency degenerate (or quasi-degenerate) case, that is routinely obtained by tuning the 
phase-matching conditions |36|, will be considered. Moreover, with no loss of generality we set Ax and Bx to be 
ordinary polarized beams and Ay, By to be extraordinary polarized [ p^ . Then the equations describing the OPO are: 

dtBx - ix[-{l + iK)Bx + lay^Bx + c'xBy + iiK^AxAy + E^] 

dtBy = 7;[-(l + ^K)By + la'yV^By + c'yB^] 

dtAx = 7x[-(l + i^x)Ax + ioLxS'^Ax + CxAy + iK^AlBx] 

dtAy ^ jy[-{l + i\)Ay + iayV^Ay + CyAx + IKoAlBx] (1) 

The coefficients ^x,y,l'x,y (effective cavity decay rates), A3,_y,A' (effective cavity mode detunings) and ax,y,a'x^y 
(diffraction coefficients) are defined in the appendix (see eqs. (|26|)). Some general remarks are worth to be made 
to show differences and common features between these coefficients and those previously defined for a "perfect" 
cavity (see for example [^Ijl^ ) . Due to the birefringence of the nonlinear crystal and the dichroism of the cavity, the 
coefficients of equivalent terms, in different equations, are slightly different, even for frequency degenerate fields. In 
fact they all depend on the relative refractive index and mirror transmittivity that are polarization dependent (see 
eqs. ( |2^ ) for details). The cavity birefringence can also cause the effective detuning coefficients, A^:, Aj, (Ai , AJ^), 
to be different, even at frequency degeneracy. Other parameters are the nonlinearity Kq (defined by eq. (pTj)) and 
the injected pump £^0 that is taken as a real number. This gives no loss of generality because it is equivalent to fix 
a common phase reference for all fields. For the sake of simplicity, we take the pump to be linearly polarized in a 
direction parallel to the phase-matched component of the intra cavity field Bx- Hence, the highly mismatch component 
By neither is pumped nor is nonlinearly coupled with other components. It is therefore very weakly involved in the 
dynamics, in spite of the linear coupling with Bx- 

The linear coupling coefficients Cx,y, c'x y account for the dichroism and the birefringence of the cavity. They read: 

p + iS / \ I p' + i6' 

Cx,y = TfTX 20 , = , ——sin{2(j)) (2 

1 ±pcos[Z(p) 1 ' ± p'cos(Zip) 

where the plus (minus) sign applies for the x (y) polarized component. Although their derivation and the exact 
relation with the physical parameters describing the cavity can be found in the appendix it is useful to describe them 
briefly at this stage. The mirror dichroism is represented by the ratio between the difference of the reflectivity and 
the average reflectivity at a certain frequency: respectively 2p for Ax^y (signal and idler) and 2p' for Bx^y (pump 
components). The mirror birefringence causes a different phase change: 2S for Ax^y and 26' for Bx^y Finally T (T') is 
the average transmittivity of the signal/idler (pump components) and is the relative angle of rotation between the 
crystal and cavity birefringence axes (the axes of dichroism are supposed to coincide with those of cavity birefringence 
for the sake of simplicity). We note that similar linear coupling terms between signal and idler, were previously 
considered by considering the insertion of wave-plates (such as quarter-wave or half-wave) in the cavity of a 

type-II OPO. In these cases the general relation Cx = — c* is satisfied. 
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III. THRESHOLD ANALYSIS AND HOMOGENEOUS PHASE LOCKED SOLUTIONS 



A. Threshold analysis 

The trivial solution of eqs. (|^) corresponds to the case in which a type-II OPO is below the threshold of signal 
generation. It is given by: 

Ax — Ay — 
Bx = cEo 

c'cEq 

TTTa; 

where c = c'' + ic* = (1 + i/^y)/[l - A^AJ^ - c'^Cy + i(A^ + AJ^)]. The threshold for instabihty is determined by 
linearizing eqs.([^) around this solution and looking for values of the bifurcation parameter Eq (the pump amplitude) 
for which perturbations grow. The general type of perturbation is given, as usual, by plane waves exp(i(f • r — Xt), 
where X{q} is the growth rate of the perturbation and q is its transverse wave vector. We first recall the main results 
of the analysis for Cx.y = c'^ y — ||l8|. In that case the trivial solution is stable for Eq < Ec, where: 



= V 1 + (4) 

and where the effective detuning A is defined as: 

~ ^ 7xA^ + 7^A^ 

Ix+ly 

If the pump amplitude exceeds Ec, the steady state becomes unstable and the signal and idler fields are generated. In 
particular, for negative effective detuning, pattern formation occurs, as studied in rcf. [^ : for this reason hereafter 
only the case A > will be considered. In this case there is a Hopf bifurcation in which homogeneous perturbations 
q — Q with 

Jm(A) =uj = ^^^^(A^ - A„) (6) 

have the largest growth rate. At threshold, a family of homogeneous oscillating solutions bifurcate from the trivial 
steady-state. For Cx^y = 0, eqs. are invariant under the transformation Ax Ax Bxp (iff)) , Ay Ay cxp {^icf)), 
and therefore the relative phase between the homogeneous oscillating solutions Ax, Ay is arbitrary. 

The introduction of the polarization coupling Cx.y ^ breaks the invariance under changes of the relative phase. 
One expects that such coupling should be able to produce phase- locked homogeneous stationary (i.e., zero frequency) 
solutions above threshold. The linear stability analysis of (^ is rather cumbersome when Cx,y ^ 0, and simple 
analytical expressions for the threshold analysis are not found. It is in this case more convenient to determine a 
threshold through the condition of existence of the relative phase locked solutions. We consider the special case 
Cx = Cy = Cx + iCx ({Cx,Cx) G TZ) which through eqs. (||) can be seen to correspond either to set cj) — ±7r/4 or 
(j) = ±3tt/4 (i.e. to fix the angle to the value that maximizes the couphng strength) or to set p = = (i.e. no 
dichroism) . We note that setting p = reproduces a particular case of the polarization coupling considered in ||l6|,|l^ . 

Let us call Ax = Ux exp{i4>x) and Ay = ay exp^iipy) the homogeneous stationary solutions; by substituting such 
formulas into eqs. |l]one gets: 

. , , , N A:, - A^ 

sin{tpy - %px 



V(c^ + 4A,)(c- -l-c^A^) 



cos(V, + i^x) = J [(A. -f T^jY - (1 + Vy - 2N/f cos(V, - ^.)(4c'- + eld)] 
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]-—;[2KoEo\/f{c'' sin{%ljy -f ijx) - c'cos{i^y + ^Jx)) + 2Vf cos{i^y - i^x) - (1 + T)] 



V 

where T = « + 4A,)/(4 -f c^A^,) 



al = Tal (7) 
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^From eqs. (|7|) we find two conditions for the existence of the phase-locked homogeneous stationary sohitions. The 
first one comes from the fact that the phase difference among the solutions (see the first of eqs. (0)) is real only if the 
modulus of the right hand side is less than one, that is if 

(A,~A,)2<4(4 + 4A,)(c; + 4A,) (8) 

This boundary in the complex plane c^,c^ defines the limit of the locking regime. Inside the boundary stationary 
solutions do not exist. Physically speaking, the locking condition means that the stationary phase-locked solutions 
exist whenever the direct polarization coupling breaking is large compared with the difference in detunings. When 
the condition (^) is not satisfied, numerical solutions show that there are still homogeneous states but their phase 
varies periodically with time. Such solutions indicate the persistence of the Hopf bifurcation found for c^^y = when 
the polarization coupling is small. 

The second condition for the existence of solutions refers to the pump value above which there is signal and idler 
generation. The threshold Ec can be determined by setting ax — into eqs. and solving for Eq. The final result 
is: 

e! = ^^2r|g|2 {(l + + (A. + TAyf + 4T\cxfcos\i:y - - 4Vfcos{^y - + T) + 4(A, + TA,)]} 

(9) 

The classification of the solutions found above this threshold is easier to understand considering the case A^: = Ay 
for which the condition (||) is automatically satisfied. In this case and regardless of the value of Cx.y the relative phase 
shift between phase-locked signal and idler can be either (in-phase solution) or tt (out-of-phase solution) and the 
amplitude of the fields is equal (ax — ay) since F = 1. Once the phase shift tpy — ipx is known it can be substituted 
into the second of eqs. (|^) that can be solved for ^y + ipx- In principle two solutions exist for each phase difference (in- 
and out-of-phase cases), due to the fact that the arcos is a multi-valued function in the range [—tt, tt]. However, if the 
negative solutions for ^y + -tpx are replaced into the third of eqs. the result is < Vi?o- Therefore only positive 
solutions of the angle sum are to be taken to guarantee that > above a certain threshold Ec. By substituting 
the phase difference and sum into the third equation the amplitude is finally found. The sign of ax can be either 
positive or negative, i.e. there are two equivalent possible solutions for the in-phase case and two for the out-of-phase 
case. The existence of these two equivalent solutions is a consequence of the symmetry (Ax, Ay) — s- ~{Ax, Ay) of eqs. 

which is preserved for Cx,y ^ 0. In summary, we find in-phase solutions (Ax = Ay) and out-of-phase solutions 
{Ax = — Ay). Each of the in- and out-of-phase cases include two equivalent solutions which we denote by the + and 
— solutions satisfying A'^y = —Ax^y. 

In the general case with A^, 7^ Ay solutions are no longer strictly in- or out-of-phase. Nonetheless, well within the 
phase-locked regime where the detuning coefficients are small compared with the strength of the polarization coupling, 
the solutions one finds are close to being in- or out-of-phase. Therefore we will still use in this situation the names of 
in- and out-of-phase solutions even if this is not generally rigorous. Other situations can occur close to the limit of 
the locking regime fixed by eq. (|^). For example, for Cx — Cy and a purely dichroic mirror, c!^ = so that F = 1 and 
signal and idler have the same amplitude ax = ay. At the onset of the locking regime | A^, — A^ |= 2 | |, and it 
follows from the first of eqs. (|^) that the two fields are locked at a phase difference ijjy — ipx = ±7r/2. In any case, for 
each locked value of the phase difference there are two equivalent -|- and — solutions satisfying A^ y = —Axy 
A bifurcation diagram for the homogeneous solutions for a generic case with A^, ^ Ay is presented in figure |l|. The 
selection of the solution that actually bifurcates, either in-phase or out-of-phase, is determined by the relative value of 
the threshold Ec for each solution. Eq. shows that the threshold is lower if the term proportional to cos{ipy — ^x) 
is positive; for the in-phase solution this occurs if > 0,c!p > and vice versa for the out-of-phase solution if 
dx < 0, < 0. An example of the thresholds calculated for the in-phase and out-ofphase solutions, as a function 
of dx — "^x: ^-^d fo'^ the same detuning values than in figure |^ are shown in figure Note that there is range of 
values of the coupling Cx,y in which eq. (|8|) is not satisfied and there is no threshold for the emergence of the phase 
locked solutions. But in the range in which (^) holds, the lower threshold decreases as the coupling strength |ca;_j,| 
increases, allowing parametric down-conversion for lower values of the external pump with respect to the reference 
case {cx,y = 0). Note also that the roles of the two solutions are exchanged if the sign of Cx,y is reversed. We finally 
note that our numerical integrations show that when switching-on the pump to a value for which both the in-phase 
and out-of-phase solutions exist, the solution which is always selected is the one with lowest threshold, while the other 
is unstable. Therefore, in practice we only find the two equivalent solutions of lowest threshold. 
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B. Polarization properties of the phase locked solutions 



An important question is the polarization state of the phase-locked homogeneous stationary solutions that we have 
just described. It is useful to consider the polarization representation given by the normalized Stokes parameters, 
defined as l33l: 



Sl 



S2 



\A.?~\Ay\'' 
\A.? + \Ay? 

Ax Ay + A*^Ay 
"12" 



v\ 



\AxV + \A 
i{AxAy — A^Ay 



*3 - ^n^^T^n^ (10) 



IAJ2 + IA 



yi 



These real parameters are sufficient to characterize any state of polarization of a monochromatic field by assigning to 
the field a point in the Poincare sphere. The equator of the sphere (53 — 0) corresponds to linearly polarized states 
and the poles of the sphere (si = S2 = 0, S3 = ±1) correspond to states of opposite circular polarization. By replacing 
the homogeneous solutions in eqs. (]To|), the state of polarization of our phase locked solutions is represented by: 

1 -r 

S3 = -J-pp'S^'^lV'y - -00;) (11) 

These equations show that the polarization state of the optical field is determined by locked value of the phase 
difference of signal and idler. For — Aj,, F = 1 and sm{ipy — ijjx) = so that (si, S2, S3) = (0, ±1, 0), where the 
plus (minus) sign applies for the in- (out-of) phase solution. This means that the two possible phase-locked solutions 
are actually linear and orthogonal polarizations whose azimuth angles arc = atan{s2/si)/2 — ±7r/4. We mentioned 
before that there are two equivalent solutions for each the in- and out-of-phase solutions. These correspond to linearly 
polarized states along the same direction, but in opposite senses, and they have the same Stokes parameters. The 
Stokes parameters are determined by the relative phase, while the two equivalent solutions have a different global 
phase. For example, in the two equivalent in-phase solutions S2 = 1 and Re(Ax) = Re{Ay), but in one of the solutions 
Re{Ax) > and in the other one Re{Ax) < 0. 

When detunings are different, the homogeneous solutions become elliptically polarized beams (33 7^ 0). However 
if detunings are small (c^A^j j, << c^) still F ~ 1 and well within the locking regime sin{'4'y — 4'x) — 0. In these 
circumstances the state of the beam is close to be linearly polarized (53 ~ 0) with azimuth angles close to ±7r/4. 
However, it is important to note that changing the detuning parameters and the strength of the polarization coupling 



it is possible to explore arbitrary states of polarization. These states are determined from eq. (11) in terms of the 
phase difference of the locked state. For example, in the case mentioned above of a purely dichroic mirror, c\. = 0, 
we have F = 1 and therefore si = 0. In this case, and at the onset of the locking regime, S2 = 0, S3 = ±1 so that 
the locked solution is circularly polarized. Going into the locking regime the polarization state will evolve towards 
linearly polarized states but keeping si = 0. 

The threshold decrease due to the polarization coupling discussed above has now a simple physical interpretation. 
Let us recall that the condition = Cy means that the relative angle of inclination (j) must be one of these values: 
[±7r/4, ±37r/4]. In particular, through eqs. (^), (c^) is positive for p > (5 > 0) and (j) = 7r/4, — 37r/4 or p < 
{6 < 0) and (p — — 7r/4, 37r/4. When A^ and Ay are in phase, the total field is linearly polarized at tt/A (or 37r/4) 
radiants with respect to the crystal axes. The conclusion that can be taken is that the intra cavity field is actually 
oriented along one of the principal axes of the cavity birefringence-dichroism. The same occurs for the out-of-phase 
solution that is a beam linearly polarized at an angle — 7r/4 (or 37r/4) radiants and then the role of coefficients is 
exchanged. So, the polarization selected is always aligned with one of the cavity principal axes. 



IV. PHASE POLARIZATION DOMAIN WALLS IN ONE TRANSVERSE DIMENSION: BLOCH-ISING 

TRANSITION 

In the previous section we have discussed the existence of two equivalent homogeneous solutions which we named 
as -|- and - solution. These are the solutions with lowest threshold, while we mentioned that other solutions of 
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higher threshold are seen to be unstable. When the OPO switches-on after setting the pump to a value above its 
threshold value, given by (^), either the + or - solution can be selected since they have the same growth rate. When 
taking into account the transverse spatial degrees of freedom, this selection, or spontaneous symmetry breaking of the 
homogeneous solution, can be local, with a different solution emerging in different spatial regions. It is then expected 
to find domain walls that separate the spatial domains with different but equivalent solutions. For either the signal 
or idler what changes when going from one solution to the other is just a sign. For example takes values A and 
—A at opposite sides of the wall. Therefore the walls can be considered as phase walls of a complex field like the ones 
described for type-I degenerate OPO When considering signal and idler the domain wall separates two 

solutions with polarization properties determined by the locking of the relative phase of the two fields. One might 
then talk about polarization walls. However, we already mentioned that the + and - solution have the same Stokes 
parameters, but there is a change in the global phase of the polarization state. In this sense we refer to these walls as 
phase polarization walls. In any case, the polarization state might present interesting features in the core of the wall. 

Phase domain walls for a complex field can be of Ising or Bloch type [|lO|,^ . As a general characterization, in an 
Ising wall there is a single field profile connecting one homogeneous solution with a second equivalent one, while we 
talk of Bloch type wall when there are two different field profiles (walls) connecting the two solutions. A Bloch wall 
implies, therefore, spontaneous symmetry breaking for the domain wall. In the following we study Ising and Bloch 
domain walls, the transition between them, and the polarization properties in one transverse spatial dimension. The 
characterization of some properties of these walls is much more clear in one dimension. Other features associated 
with two-dimensional phenomena are postponed to the following sections. 

Numerical integrations | ]37[ | of eqs. (|l|) confirm that stationary uniform domains of the + or - solutions form 
spontaneously starting from a randomly and weakly perturbed trivial steady-state (|^). Well within the locking 
regime, the domain walls are of the Ising type, but changing the values of c^.y and moving towards the boundary of 
the locking regime we find a transition from from Ising to Bloch domain walls llQ] . An example of an IW, is presented 
in figure It connects the -I- solution at x ~* — oo with the -I- solution at a; — > oo. By plotting the numerically 
obtained solution in the complex plane {Re{Ax), Im{Ax)) (figure ||b) it is clear to observe that the IW is characterized 
by a zero crossing of the field. An example of ID optical BW is instead given in figure U, a) and b); note that the field 
amplitude (represented by the vector modulus in the complex plane) never goes to zero and the wall consists of an 
almost pure phase rotation of tt radiants. The phase can rotate in two possible senses along the interface, clockwise 
or counterclockwise in the complex plane. This characteristic is usually called the wall chirality and it is defined to 
be positive for clockwise rotation, or negative for counterclockwise rotation. Therefore, there exist two equivalent 
domain walls of opposite chirality for A^- One of the two appears by spontaneous symmetry breaking. In the example 
of figure ^ the wall of negative chirality for A^ is selected. 

In the example of the Ising wall of Fig. |3| parameters are such that we are well within the locking regime and the 
domain wall connects homogeneous equivalent in-phase solutions in which the locked phase difference is close to zero 
{ipy — ipx — —0.122). The shape of the field Ay across the domain wall is similar to the one of A^- The situation is 
different for the example of Bloch wall in Fig. ^. Parameters correspond here to a situation close to the boundary 
of phase- locking and the homogeneous 4- and - solutions have a locked phase difference ipy — ipx = it /2. In addition, 
Bloch walls for this system are characterized by the fact that the wall profile for the field Ay has always opposite 
chirality to that of A^, as seen in the example of Fig. ^b. 

The polarization characteristics of Ising and Bloch domain walls are very different. For an Ising wall the Stokes 
parameters are seen to remain constant across the core of the wall. This is due to the fact that the phase difference 
i'y — i^x remains fixed to its locked value while going from the - to the -I- solution across the wall. On the contrary, in 
a Bloch wall the locked value of the phase difference ij^y — V's is a function of the position while moving form one side 
to the other of the wall. The consequence is that the Stokes parameters have a nontrivial space dependence across 
the wall determining peculiar polarization characteristics of the core of the wall. As examples of such polarization 
characteristics we show in Fig. |^ the variation of the Stokes parameters across two examples of Bloch domain walls. 
In Fig. I a we consider a Bloch wall which connects two linearly polarized states. In Fig. ||b, which corresponds to 
the wall of Fig. |[ the wall connects two elliptically polarized states which are close to being circularly polarized. 
In Fig. ^ the two asymptotic states for x ±cxd are in-phase solutions characterized by (51,52,^3) = (0,1,0) 
corresponding to a linearly polarized state of azimuth 9 — 7r/4. At the core of the wall (si, S2, S3) — (0, —1, 0) which 
corresponds to a state of orthogonal linear polarization. Note that as the phase solution rotates tt radiants the Stokes 
parameters si, S2, S3 return to the initial values, i.e. the polarization of the total field is the same on each side of the 
wall, as previously mentioned. Along the wall the polarization changes, the field becoming elliptically polarized, but 
not in an arbitrary manner. The transformation is forced to occur for Si ~ 0, i.e. the azimuth angle of the ellipse 
is practically fixed to 6 ~ ±7r/4. For cc < 0, and close to the core of the wall, S3 ~ —1, indicating a state close to 
being left circularly polarized, while for a: > 0, and close to the core of the wall, S3 ~ 1, indicating a state close to 
being right circularly polarized. In terms of the Poincare sphere, the representative point moves from a point in the 
equator through the vicinity of the south pole to the opposite point in the equator, and back to the original point 
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through the vicinity of the north pole. The change of eUipticity across the wall, ri{x) = arcsin(s3)/2, yields a natural 
interpretation of the chirality: a BW of negative chirality, like the one shown here, means that the eUipticity of the 
polarization state first decreases to a maximum negative value going as we move to the other side of the wall to a 
maximum positive value. For a BW of positive chirality the excursion in eUipticity goes in the opposite direction. 

In Fig. the variation of the Stokes parameters indicates a sequence of cUiptically polarized states with two points 
in the core of the wall at which the state becomes linearly polarized (53 = 0). The change in polarization state still 
occurs for si ~ because being = 0, still F = 1. The representative point moves now in the Poincare sphere from 
a point close to the north pole to the vicinity of the south pole crossing the equator, and back to the original point 
along the other side of the sphere crossing again the equator. An opposite sense of motion along the sphere would 
correspond to an opposite chirality of the wall. 

We note that a quantitative precise description of the variation of the Stokes parameters can be generally given by 
invoking the relation, that follows from symmetry considerations, ~ lA* point wise along the wall. The Stokes 
parameters as functions of the spatial variation of the phase ipx{x), are then given by (substitute A^ = iA* into eqs. 
®: 

si{x) ~ 

S2{x) ~ sin{2tljx{x)) 

S3{x) ~ -cos{2Mx)) (12) 

We next turn to consider the dynamics of the domain walls, which is also useful to determine the transition between 
Ising and Bloch walls. Isolated Ising walls in ID are stable and remain stationary. The dynamics of ID BWs depends 
critically on the values of the cavity decay rates and the detuning. For '^x^x = ly^y^ ID BWs do not move: they are 
stable stationary interfaces between equivalent uniform domains. A similar situation takes place in the potential limit 
of the ID parametrically forced Complex Ginzburg Landau Equation (PCGLE) where stationary BWs have been 



found analytically |21 3q|. On the contrary, for ^x^x 7^ ly^y walls of different chirality move in opposite directions, 
as it also happens outside the potential limit of the (PCGLE). The velocity of the resulting BW depends on the value 
of the parameters, in particular it depends strongly on Cx,y In figure || the velocity of BWs as a function of Cx (real), 
for selected values of the other parameters, is shown as it results from numerical solutions. For small Cx BW are not 
stable because we are outside the locking regime ( eq. (^ is not satisfied); for larger values the velocity decreases 
by increasing Cx and finally it vanishes. The vanishing of the velocity identifies the transition point in which a BW 
decays into an IW. This transition is continuous, i.e. the amplitude of the wall solutions, that is almost constant for 
BWs for small Cx, shows larger and larger variations as the transition is approached. The amplitude of the signal 
field (and also the idler) become small at the core of the wall and it eventually reaches a zero. At this point BWs and 
IWs exchange their stability and only standing IWs are found beyond that critical value. The conclusion is that in 
the regime of the parameter space for which BWs are stable, IWs are unstable and vice versa. 

In addition to the transition from Bloch to Ising type, the strength of the direct polarization coupling coefficients 
Cx,y also controls the wall width. Moving well within the locking regime the width becomes small, while it diverges 

fiS Cx^y ^ 0. 

We finally mention that other forms of stable localized structures or domain walls are sustained in this system, although 
they do not appear spontaneously when starting from a weakly perturbed trivial steady-state (^) and a pump value 
above threshold. They can be formed far above threshold mediated by the unstable phase locked solutions of higher 
threshold discussed in Sect. III. A first example, shown in Fig. is a stable ID localized structure in the background 
of a + in-phase solution. This is structure is generated from a step-like initial condition, shown in Fig. which 
connects the stable -I- in-phase solution with the unstable -I- out-of-phase solution. The dynamical evolution of this 
unstable configuration leads to the localized structure. A second example shown in Fig. [7|c is a new kind of domain 
wall which connects, as in the Ising and Bloch walls discussed above, the equivalent -I- and - solution of the phase 
locked solutions of lowest threshold. The difference with our previous Ising walls consists in a more complicated 
structure of the core of the wall. In this case the field Ax vanishes at three points in the core of the wall. This new 
type of domain wall emerges dynamically from a similar initial condition that the one shown in Fig. |^b, but changing 
the sign of the unstable solution of the step: The initial step like condition connects the stable -|- in-phase solution 
with the unstable - out-of-phase solution. 



V. CHARACTERIZATION OF 2D BLOCH WALLS 

In our 2D numerical solutions [ ^Tj with random initial conditions around the trivial unstable solution we also observe 
the spontaneous formation of BWs which are now lines in the transverse plane. These walls evolve dynamically as 
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described below. A main new feature of BWs in 2D is that the domain walls can emerge with an opposite chirality 
in different spatial regions along the wall. The points on the wall where the chirality changes sign are singular points: 
At these points the phases of the signal and idler field are not defined and the amplitudes go to zero, i.e. they can be 
classified as topological defects. The BW, in that particular point, actually degenerates into an IW. A snapshot of 
typical transient transverse patterns generated by the Bloch walls is shown in figure |8| for the component A^. The two 
equivalent + and - phase locked homogeneous solutions are represented in Fig. ^ a) by regions of different intensity on 
a gray scale (A^). Likewise, the segments along the walls of different chirality are represented respectively by black 



or white segments {B± ) in figure 
dots in the intensity field (figure 



^ a). The defects, where the changes of chirality take place can be observed as black 
8 b)). Note that, except for the defects, the field intensity is almost constant in the 
domains and only slightly modulated close to the core of the wall. The phase field, shown in figure || c), demonstrates 
that phase defects with topological charge ±1 occur at the points of change of chirality on the wall. An expanded 
vision of the amplitude of the field in the vicinity of one of these defects is shown in Fig. ^. The corresponding 
snapshot for the field Ay shows the same structure with walls and defects at the same positions as for A^- 

There are different effects determining the dynamics of BWs in 2D. A first mechanism is the one found in ID, related 
with detuning and damping coefficients. However, in 2D the dynamics is also influenced by the curvature of the walls 
and the presence of field defects. In what follows two main regimes are identified. In a first regime, called of domain 
growth, curvature effects dominate. This leads to the complete disappearance of all fronts, all defects and all domains 
except one, the final state being an homogeneous solution. The second is characterized by a persistent dynamics in 
which defects are deeply stable objects, while domains and walls are continuously generated and annihilated. 



A. Regime of domain growth 



For jx^x = ly^y, flat 2D BWs are stable structures. This corresponds to the fact that ID walls, for the same 
values of the parameters, do not move. The transient dynamics is then ordering process mainly controlled by the 
curvature effects in which the walls evolve by reducing curvature. This leads to the growth of one of the two equivalent 
solutions at the expenses of the other and the annihilation of all the defects. This process is shown in figure where 
snapshots at different times are shown. The images of the upper row show the intensity of polarization component 
Ax , while in the pictures of the lower row, the real part of the same field is presented. Slowly all the defects and walls 
disappear and the final state is a homogeneous phase-locked solution. 

In this regime the normal velocity of the fronts is determined by the local curvature of the wall. This is demonstrated 



by the result of figure 11 where the evolution of the square radius of a domain surrounded by a circular BW has been 
determined. The result is a growth law R{t) ~ t^/^ characteristic of curvature driven domain growth [ p6[ . Similar 
structures, annihilation dynamics and dynamical exponents for the growth law of BWs domains, have been reported 
in the description of the ordering process of a non conserved anisotropic XY-spin system in 2D [p8| . 



B. Regime of persistent dynamics 

More peculiar of our system is the reginie found for ^x^x 

jyAy. Let's recall that in this regime ID optical BWs 
of different chirality move in opposite directions while IWs (characterized by a point of zero amplitude) do not move. 
The dynamics in 2D is reminiscent of this behavior; in fact the defects are notably stable as the corresponding ID 
IWs, while BWs of different chirality move in opposite directions. The combination of these two effects is that BWs 
spiral around the defects; such phenomenon was also observed for BWs in other physical systems |39|. The spiraling 
dynamics of an isolated defect is shown in Fig. |l^ obtained with a flat-top profile for the pump beam. The stability of 
the defects, with these physical boundary conditions is remarkable. This is an important result because stabilization 
of optical vortices has been always critical in lasers and nonlinear optical systems pO|] . When many of these defects 
arise spontaneously from random initial conditions around the trivial unstable solution, the system becomes trapped 
in a complicated state of persistent dynamics in which an homogeneous state is never reached. An understanding 
of such persistent dynamics is easier for particular initial conditions. In the example shown in Fig. (|l^) the initial 
condition is a flat BW with equally spaced defects along the wall. This gives rise to segments of different chirality 
along the wall that start moving in opposite directions. The net result, after collision of spiraling BWs of the same 
chirality is the periodic emission, in each direction, of BWs of alternating chirality. The defects remain stable and 
walls are re-generated by the spiraling process. The dynamical process that we have just described becomes fuzzy 
close to the Ising-Bloch transition. In such regime of parameters the amplitude of the field becomes very small at the 
core of the wall and the clear distinction between point defects and wall segments of different chirality is lost. 



9 



VI. 2D ISING WALLS 



Beyond the transition from BWs to IWs the latter appear spontaneously separating spatial domains of the + and 
- solutions, but their dynamics and the asymptotic state of the system still depends strongly on the strength of the 
direct polarization coupling. As in the case of BWs, we find for IWs two main regimes: a regime of domain growth 
and one of labyrinthine pattern formation. 

For coupling values close to the Bloch-Ising transition, a flat Ising wall is stable. The transient dynamics is then 
controlled by the curvature of the walls. Curvature reduction leads to domain growth much in the same way hat we 
already described for BWs in the corresponding regime. In figure |4| a typical time series of snapshots of the transverse 
Ax field is shown. The final asymptotic state is homogeneous. Similar dynamics has been already reported in type-I 
DOPO above threshold |]8|, |ll| , ^ , the only difference being the vectorial character of the field in this case. 

Far from the Bloch-Ising transition and moving deeper into the locking regime transverse labyrinthine patterns 
are spontaneously formed in the system starting from a random perturbation of the trivial unstable solution |p6| , p7| . 
Snapshots of the evolution of a pattern of this type is shown in intensity and real part in figure |l^. Note that the time 
evolution of labyrinthine patterns is very slow. The creation of a labyrinthine pattern stems from the fact that, in 
this regime of parameters, flat Ising walls (i.e. with no curvature) are modulationally unstable. Roughly speaking the 
finger growth is associated with a band of modulational frequencies of the front curvature that tend to increase their 
curvature. This is reminiscent of what has been reported for intra cavity second harmonic generation and vectorial 
Kerr resonators pGf . In order to illustrate the modulational instability, the evolution of an initially perturbed Ising, 
flat wall is shown in figure |l^. 



VII. CONCLUSION 



In conclusion we have demonstrated that Bloch walls can be found in nonlinear optical systems, in particular in 
type-II optical parametric oscillators. They appear when there exists a small birefringence and/or dichroism of the 
cavity along axes that do not coincide with those of the nonlinear crystal. These effects introduces a linear coupling 
between the signal and the idler which causes self phase-locking of the two fields. There exist two possible steady-state 
solutions, characterized by a phase shift of n radiants of both polarization components and thus different domains 
spontaneously form in which one or the other solution is selected. Separating walls can be either of the Bloch or the 
Ising type depending on the strength of the coupling coefficient. For small values Bloch walls are stable and appear 
spontaneously above a predicted threshold out of a random perturbation of the trivial steady-state. Bloch walls have 
been characterized both in one and two dimensions. In one dimension a physical interpretation of Bloch walls is given 
in terms of polarization variations that connect two homogeneous states that represent the same state of polarization. 
The chirality is instead related to the ellipticity variations. Multiple hump Ising walls have also been found starting 
from particular initial condition. In two dimension, Bloch wall peculiarity is the possibility of having wall sections 
of different chirality, i.e. where the phase rotates in two possible ways, clockwise or counterclockwise in the complex 
plane. Where chirality changes sign the phase has defects, where the field amplitude is zero, and the wall degenerate 
into an Ising one. Two dynamical regimes, that depend on the decay rates and the detunings, are found: in the 
first one, the wall dynamics is dominated by the curvature and a final homogeneous state is reached; in a second 
regime, the walls spiral around stable defects and a persistent creation and annihilation of fronts is observed. The 
transition from Bloch to Ising walls has been also observed when the linear coupling strength is increased; transition 
is characterized by larger and larger variations of the amplitude of the field close to the core of the wall and the 
fact that walls stop moving. Ising wall dynamics has been also considered, in particular their curvature modulation 
instability, that leads to the creation of a labyrinthine pattern. 

This work has been partially supported by the Spanish MCyT project BFM2000-1108 and by the European Com- 
mission project QSTRUCT (FMRX-CT96-0077). The authors acknowledge clarifying discussions with G-L Oppo. 



VIII. APPENDIX 

In this appendix the derivation of eqs. (|l]) is presented. The main difference with respect to previous derivations 
of mean filed equations for OPOs is the inclusion of the effects of dichroic and birefringent mirrors. In order to 
simplify the model, let us assume that only one, out of the four mirrors of the ring cavity, is birefringent and dichroic. 
This means that, in a proper system of orthogonal axes, the matrix M, that represents the transformation of two 
orthogonally polarized component of a beam, in the Jones formalism p2], is: 
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r2e-^^'- J - [ {l-p)e-^S ) (1-^) 

where r = {ri + r2)/2,p — (ri — r2)/(2r),(7 — —{ipi + V'2)/2 and 6 = ('02 — '0i)/2- The dichroism, i.e. the different 
reflectivity of the mirror for different polarizations (ri ^ r2), implies that p ^ 0, while the birefringence is represented 
by the fact that ■01 7^ V'2, i-e. the phases of the reflected components undergo a different change {5 0). 
A first important remark to make is that the mirror anisotropy axes might not coincide with the axes of the nonlinear 
crystal, that is also birefringent in order to realize the phase matching between the pump and the generated fields. 
In other words the anisotropic crystal has its own preferred polarization axes that can be rotated by an angle with 
respect to the mirror principal axes (those for which the mirror matrix is diagonal). Therefore, when passing from 
the propagation in the cavity axes reference frame to that in the crystal a rotation, represented by the matrix 

\ striq) coscp J 

has to be applied. At the end of the propagation in the nonlinear medium an inverse transformation (i?^^) is needed 
to restore reference frame of the cavity axes. 

The generic n + 1 round trip in the cavity for the signal and idler vector field of the previous step n (-E„) is then 
represented by the following transformation (all mirrors except the last one are perfectly reflective, i.e. their matrices 
are all equal to the identity matrix): 

^„+i = M /(i?^„) (15) 

where /(•) is the result of the propagation inside the nonlinear medium. A similar formula can be written also for the 
pump vector field Fn that, in principle, has both polarization components. Usually one component does not participate 
in the nonlinear dynamics (it is not phase-matched) and it is neglected a priori. In this case, due to polarization 
coupling, it is included in the model although the final result is that, under not very restrictive hypotheses, its effects 
can be neglected. 

The output of the function /(•) is the vector field as it results after the integration of a set of nonlinear ly coupled 
propagation equations, i.e. it involves products of En and Fn- Hereafter only the signal and idler vector field £"„ 
equations will be considered; similar calculations can be repeated for the pump components. 

Let's remark that the optical carrier frequency of each component has been already removed (i.e. envelope equations 
p6| are searched for) and that the carrier frequencies of all waves are determined by three conditions: phase-matching 
and energy conservation of the nonlinear interaction and the condition of resonance due to the cavity. It has been 
demonstrated that, for a type II OPO, there are several signal-idler pairs of oscillation frequencies that can satisfy 
these conditions. Among these solutions there is also the case of quasi, or totally, frequency degenerate signal and 
idler. This is also verified experimentally by the fact that type-II OPO's, differently from type-I, can be smoothly 
tuned through frequency degeneracy jl^^Q. Let us remark, finally that for type 11 OPOs signal and idler are always 
polarization non-degenerate. 

The boundary condition to impose into eq. ( [l5| ) for steady-state operation of the OPO is that the round trip 
transformation coincides with the identity, i.e. En+i — En- Let us define A{L) — f{REn) the vector field at the output 
of the crystal of length L (that is also the cavity length for the completely filled cavity) and A{0) = REn ~ REn+i the 
field at the crystal input; by multiplying the left and right hand side of eq. (|l^) by R (on the left) and substituting 
previous definitions the following equation is found: 

1(0) = RM R-^ A{L) (16) 

Let us now set A'{z) = H{z)A{z) such that 

A'{0) = 1(0) , A'iL) = 1(0) (17) 

the second condition imposing the periodicity after a round trip. The general form of the matrix H{Z) satisfying ( p^ ) 
is: 

H{Z) = ( ^hy^^ _ ^-hy^Z ^hyyZ ] (l^) 



This matrix extends the scalar transformation used by Lugiato and Oldano in their original paper |41|, devoted to the 



study of stationary spatial patterns in optical systems with two-level atoms, to a vectorial case. It is easy to verify 



11 



that H{0) is the identity matrix while H{L) = RMR ^ and thus all the elements of the matrix hij can be explicitly 
calculated: 

hxx — —{Inlr) — icT + ln[[l + pcos2(j))cosS + i{p + cos2(l))sin8W 
L 



hyy — —{ln{r) ia + ln[{l — p cos2(f))cos5 + i{p — cos2(f))sinS\} 



hxy — ~ —{ln[[pcos5 + i sin6)sin2(j) + 1/ {pcos6 + i sinSY si'n?2(j) + 4] — ln{2)} (19) 
ij 



The evolution of the field in the crystal is governed by equations of the type: 

d,A = C{A)+N(A,B) (20) 
where B = [B^, By] = RF^ represents the pump vector field. The linear term is 

and includes the diffraction {k^^y are the longitudinal wave vectors of signal and idler, is the spatial transverse 
Laplacian operator) and the phase velocity mismatch {v^^y are the phase velocity respectively of signal and idler, dt 
is the differential operator with respect to time). The nonlinear operator is 

U{A,B)^iKBx(^l l^A* (22) 

where K is the nonlinear coefficient. Since A = H^^A' the evolution of the field A' in the crystal can be determined 
from: 

d,A = {d,H-^)A' + H-^d,A' (23) 

that finally yields: 

d,A' + H{d,H-^)A' = HCH-^A' + HJ\f{H*)-^A' (24) 

The matrix H is known and so all the terms of the last equation can be explicitly calculated; in particular, by following 
the guidelines of ref. |4^, if p and 5 are small parameters (of the order of the transmittivity T=l-r of the mirror) 
one can calculate a set of first order (in T) equations by expanding all coefficients up to this order. After long but 
straightforward calculations, the final result is a set of coupled equations for the vector A' — [Ax,Ay]^ , that still 
contains both and dt operators, but includes the boundary conditions. It is worth to write these equation because 
the coefficient of each term appears explicitly written in terms of physical parameters: 

dtAx + VxdzAx = hxxVxAx + ^{p + iS)sin2(f)Ay + i^^\/'^Ax + iVxKA*Bx 

dtAy + VydyAy = hyyVyAy + ^ {p + iS) sin2(j)Ax + i^^'^Ay + iVyKA^Bx (25) 

By exploiting the single longitudinal mode approximation, which is quite a good one for continuous wave OPO's, the 
longitudinal spatial dependence can be finally removed and the equations describing the time evolution are exactly 
eqs. (|^). The coefficients of eqs. (|l|) can be found from eqs. ( |25| ) once the expressions ( |l9| ) are substituted. They 
read: 

Vx y{T ±pcos24>) 

^ _ Vx,yio' ± S COs2(t)) 

"''^ ^ T±pcos2(j) 
L 



"^'^ 2kx.y{T ±pcos2(i)) 

Cx.y = 7F^-^-^s«"-2(/) (26) 
1 ±pcos2(p 
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where the plus (minus) sign apphes for the x-polarized (y-polarized) component. Note that all coefRcients can be 
different because of the crystal birefringence (vx 7^ Vy), the cavity birefringence and/or dichroism (p 7^ 0, (5 7^ 0). 
Finally the nonlinear coefficients is defined as 

KL , , 

K^ = — (27) 

Actually it is also slightly different for the two polarization, but this difference has been neglected because it is only 
due to the mirror dichroism; all the other coefficients have larger differences because two effects (crystal and cavity 
birefringence) contribute. 

Regarding the coupling coefficients some special cases are worth to be remarked in relation to what is discussed in 
other sections of this paper. For example when cos2(j) = 0, the special condition — Cy (complex) is obtained; in this 
case (j) — ±7r/4, ±37r/4, sin2(j) = ±1 and therefore the linear coupling among polarization is maximized in modulus. 
Other two special cases are the purely birefringent mirror p = 0, for which — Cy are purely imaginary, and the 
purely dichroic mirror that yields 6 = and thus purely real Cx,y 
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FIGURES 



FIG. 1. Bifurcation diagram of the homogeneous solutions (soUd curves stable, dashed curves unstable). Parameters are 
Ko = l,Cy=c^= 0.05 (1 + i), c'^ = c'y = 0.1, A'^ = Ay = 0, = 0.01, Ay = 0.03. 

FIG. 2. Threshold of instability Ec for the trivial stationary solution (eq. (^): solid (dashed) curve is the threshold of the 
in- (out-of-) phase solution. Here Cy = Cx and the others parameters are the same of figure ^ 

FIG. 3. Numerical solution of eqs. (^) in one spatial dimension showing and example of an Ising wall: a) Solid (dotted) 
curve is the real (imaginary) part of the field Ax as a function of the transverse coordinate; b) the same wall is represented in 
the complex plane {Re{Ax), Im{Ax)) in solid line. Dotted curve is the Ising wall associated to the field Ay. The parameters 
are 7y = 7y = 1, Jx ^ I'x = 1-001, A'x = A'y = 0, A^ = 0.01, Ay = 0.03, a'x ^ a'y = 0.125, ax ^ ay = 0.25, Kq = I, Eq ^ 1.25, 
c'x = c'y= 0.01(1 + i) and Cx = Cy = 0.082. 



FIG. 4. Numerical solution of eqs.(y) in one spatial dimension, showing a Bloch wall, a) Solid (dotted) curve is the real 
(imaginary) part of the field Ax as a function of the transverse coordinate; b) the same wall is represented in the complex plane 
{Re{Ax), Im{Ax)) in solid line. Dotted curve is the Bloch wall associated to the field Ay. Parameters are the same of figure 
except Cx — Cy = 0.01. 



FIG. 5. Variation of the polarization along a Bloch wall, represented by means of the Stokes parameters (^^ as they result 
from a numerical solution: si{x) (solid curve), S2{x) (dashed), 53(2:) (dotted), a) The Bloch wall separates linearly polarized 
states. Parameters are the same of figure 0, except "fx = Jx = li A3, = Ay = 0.03, c'x ~ c'y — 0.01, Cx = Cy = 0.02; b) The same 
representation for the Bloch wall of Fig. ffwhich separates elliptically polarized states. 

FIG. 6. Velocity of a BW as a function of (c!^ = 0) as given by numerical simulations of eqs. (^ in the regime jxAx 7^ JyAy . 
The parameters are the same of figure^ except the values of Cy — Cx. For the last three points in the figure (corresponding to 
Ising walls) the velocity is zero. The dotted line indicates the onset of the locking regime. 

FIG. 7. Examples of multiple humps stable solutions: in a) the field comes back to the initial state; b) initial condition used 
in (a). Here Re{Ax) = Im(Ax) for x < 0; c) IW with two crossings of the zero amplitude point. Parameters are ax = 0.25375, 
ay = 0.24625 a'x ^ a'y ^ 0.125, A^ = A^ = 0, Ax ^ Ay = 0.02, 7^=7^, = 1, jx = 0.985, -jy = 1.015, c'x ^ c'y ^ 0.01, 
Cx = Cy = 0.2 + i 0.02. 



FIG. 8. A snapshot at time t=1600 of the field Ax{x,y,t), spontaneously generated from random initial conditions close to 
the trivial steady-state. Figure a), b) and c) show respectively the real part, the intensity and the phase. Black and white 
segments in the walls {B±) in a) correspond to opposite sense of rotation of the phase (chirality). The black points in b) are 
the defects, where signal and idler amplitude vanishes. The parameters are the same of figure]^ except A^, = 0.01, Ay = 0.03, 
c'x ^ c'y ^ 0.025(1 - i/2) and Cx = Cy ^ 0.02(1 + i). 



FIG. 9. Amplitude of the Ax field as a function of the transverse coordinates {x, y) in the vicinity of a defect. The parameters 
are the same of figure except 7a; = 7^ = 1, Cx — 0.02(1 -I- i/2), Cy — 0.02. 

FIG. 10. The domain growth regime {'yxAx ~ ^yAy) is presented by means of snapshots at different times: a,)t — 200; b) 
t = 600; c) t = 1000; d) t = 2000. The upper row shows the evolution of the intensity of Ax{x,y,t) while in the lower row 
the real part of the same field is shown. The initial condition is random and the parameters are the same of figure ^, except 
Ay = 0.01002, 71 = 1.002, Cx ^ Cy = c'x = c'y ^ 0.02(1 + i). 
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FIG. 11. The asterisks shows the time evolution of the square radius of a circular BW domain, in the growth domain 
regime, as they result form the numerical solution. The solid line is a linear interpolation of numerical data. Parameters are 
the same of figure |l^, except c'^ = c'y — 0.01(1 + i) and Cx — Cy — 0.02. 

FIG. 12. Time evolution of one spiral of BWs around a defect. In the left column the amplitude of Ax{x, y, t) is represented 
while in the right column the real part of the same field is shown, a) t = 560; h) t = 1120. The parameters are the same of 
figure |. 

FIG. 13. Time evolution of a BWs, spiraling around an array of defects, imposed as an initial condition: a) t=0; b) t=1000; 
c) t=1150; d) t=1550, in the left column the amplitude of Ax{x, y, t) is represented while in the right column the real part of the 
same field is shown. The parameters are the same of figure |^ except A.x = 0.01, Ay — 0.03, Cx = Cy — 0.025 and c'^ = c'y — 0.01. 

FIG. 14. Snapshot at different times, observed in the coarsening regime: a.) t = 300; h) t = 900; c) t = 1900; A) t = 3000. 
The upper row presents the intensity pattern and the lower row Re{Ax). The initial condition is random and the values of the 
parameters are the same of figure ^ except: Cx = Cy = 0.09. 

FIG. 15. Formation of a labyrinthine pattern in the regime for which IWs are stable but their curvature is modulationally 
unstable: the upper row shows the intensity pattern and the lower row Re{Ax) for a) t = 75; h) t = 1500. The initial condition 
is random and the values of the parameters are the same of figure |l^, except: Aj, = 0.01001, Cx = Cy = = c'y = 0.3 i. 

FIG. 16. Modulational instability of an initially fiat IW: upper row shows \Ax{x,y,t)'i^ and the lower row the real part of 
the same field. In this regime ID Ising front are stable. The values of the parameters are the same of figure hsi. 
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